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Abstract 



We study the properties of dijet production in deep inelastic scattering using a unified 
BFKL/DGLAP framework, which includes important subleading ln(l/x) contributions. 
We calculate the azimuthal decorrelation between the jets. We compute the cross section 
for dijet production as a function of and the jet transverse momentum, as well as 
calculate the total dijet rate. We compare the predictions with HERA data. 



Dijet production in high energy deep inelastic electron-proton scattering can expose proper- 
ties of small X behaviour in QCD, as can be seen from Fig. 1. In the dominant 7* (7 qq dijet 
subprocess, the incoming gluon can have sizeable transverse momentum accumulated from dif- 
fusion in kT along the gluon chain ^. The value of the transverse momentum, and hence 
the azimuthal decorrelation between the jets, increases with decreasing x. That is the jets are 
no longer back-to-back since they must balance the appreciable transverse momentum kr of 
the incoming virtual gluon. The azimuthal decorrelation from the back-to-back configuration 
= vr is therefore a measure of kr and may be expected to be an indicator of the diffusion along 
the BFKL chain. Clearly to obtain a reliable measure we must avoid the infrared region kx — 
(that is ~ tt). However, as will be seen, we are able to make an essentially parameter- free 
prediction of the integrated dijet production rate (at the parton level). 

The description of dijet production is based on the unfolded kx factorization formula for 
structure functions [0, ^, which exposes the unintegrated gluon distribution fi^Xg^k"^) more 
locally than the structure functions themselves. The calculation goes beyond the conventional 
(fixed order) QCD-improved parton model, which is known to underestimate the dijet rate 
0. We will find that the approach driven by the unintegrated gluon distribution gives an 
enhancement in the predicted rate and hence an improvement in comparison with the data. 

The differential cross section for producing two jets of transverse momenta pir and p2T in 
deep inelastic scattering is 



da 47ra^ 



dxdQ'^dcpdplrpdplx xQ"^ 




dFr . dFr, 



d(f)dp\j'dp2'j' d(j)dp{j,dp2T 



(1) 



where, as usual, the deep inelastic variables = — g^, x = Q'^/2p-q and y = p-q/p-pe where pe,p 
and q are the four momenta of the incident electron, proton and virtual photon respectively, see 
Fig. 1. The differential structure functions are obtained from the kx factorization prescription 
by unfolding the integrations over p\j,, p\rp and the azimuthal angle between the q and q jets. 

It is convenient to use a Sudakov decomposition of the jet four momenta 

Pi = (1 - P)q' + aiP + Pit 

(2) 

P2 = Pq' + a2P + P2T 
where q' = q + xp and p are the basic lightlike momenta. Since the jets are on-mass-shell 

where is the mass of the quark. The kx factorization formula for the differential structure 
functions is 



d(j)dplrpdp2x 



I d(3 pIt, pIt, 0, Q') (4) 

q Jo rCj^ 



1 



with i = T,L. The function f{xg, k^) is the unintegrated gluon distribution describing the 
gluon chain in Fig. 1. The variables Xg and /c^ are the longitudinal momentum fraction and 
the transverse momentum, relative to the proton, carried by the gluon which couples to the qq 
jet pair. They are given by 

Xg = X + ai + a2 ( ^ [1 + 4:PiT/Q'^]x) 

(5) 

kj, = + + 2piTP2TCOS(f), 

see Fig. 1. The functions J-'i, which describe the virtual photon-gluon fusion subprocess 'jg — »• 
qq, are 



= e as{k 



8n 



Pit , P2T , '2piTP2T cos ^( I 1 , . 



where Cq is the charge of the quark and where the denominators 

A = pIt + + /5(1-/3)Q2. (8) 

The unintegrated gluon distribution f{xg,k'^) is obtained by using a unified BFKL/DGLAP 
formalism, which includes sub-leading \ia{l/xg) contributions^, to fit to the inclusive F2 structure 
function data In this way f{xg, k^) is determined for k^ > k^ where /cq = 1 GeV^. Formally 
the integration limits for (3 in (§) are and 1, but they are constrained by the condition Xg < 1. 
Moreover the lower limit on k^ in the determination of the unintegrated gluon means that we 
cannot predict the azimuthal decorrelation between the jets in the near back-to-back domain 
0^7?. Our decorrelation predictions are limited to the region 

1 + COS0 > kl/2pl (9) 

where Pq is the minimal value of the transverse momentum of an outgoing jet in the dijet 
system. 

Predictions for the azimuthal decorrelation are shown in Fig. 2. They are compared with 
the measurements made using the ZEUS detector which so far have been presented in the thesis 
of Przybycien We use the parton level data obtained with the kx jet-finding algorithm. 
The predictions use the same cuts as the data; that is > 8 GeV^, outgoing electron energy 
< 10 GeV, pr(jet) < 4 GeV, -2 < ■r]{\ab) < 2.2 and ijCRCM) > 0, where the pseudorapidities 

-"^Next-to- leading order (NLO) ln(l/x) corrections are known to be large 0, and imply the need for 
an all-order resummation. We are able to include a major part of the all-order resummation by imposing a 
consistency constraint on the BFKL kernel Q. The constraint requires the virtuality of the exchanged gluons 
along the chain to be dominated by their transverse momentum squared. Although the constraint contributes 
at all orders, its NLO contribution gives about 70% of the exact NLO result for the BFKL Pomeron intercept. 
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7] refer to the HERA and hadron centre-of-mass frames respectively. In each x bin the data for 
the distribution are normahzed to unity, and the predictions are normahzed to the data point 
centred at = 2.55 radians (0 = 146°). We see that there is satisfactory agreement between 
the predictions and the shape of the observed distributions. It should, however, be noted 
that, due to the cuts, eq. (||) implies that Xg > 9x. Thus the observed x bins do not fully expose 
the small Xg domain. Hence the broadening of the azimuthal decorrelation, although visible in 
the predictions, is quantitatively marginal. The Xg distribution of the data is compared with 
the predictions in Fig. 3(a). We see these data sample the gluon in the region Xg ~ 10^^. Not 
surprisingly, in this Xg domain the distribution does not give a definitive test of the underlying 
dynamics. Indeed Monte Carlos, which do not embody BFKL effects, can also describe the 
distribution reasonably well [Q. 

We emphasize that the calculation of dijet production is essentially parameter free. More- 
over it applies to the full kinematic domain since it is based on an unintegrated gluon which is 
obtained from a unified BFKL/DGLAP approach. The parameters which enter the determina- 
tion of the unintegrated gluon are completely specified by the fit to the F2 structure function 
data Thus we can make an absolute comparison with the measured cross section for dijet 
production. Fig. 3(b) shows the comparison as a function of Q^. At large there is excellent 
agreement. However as decreases the prediction, with its weaker dependence, falls below 
the data. The reason is that for <^ Apfj. the denominators Di of (||) are dominated by pfj, 
and hence the calculated cross section depends only weakly on Q^. There is a natural expla- 
nation of the discrepancy in Fig. 3(b). Dijets may also be produced in the photon hemisphere 
from the higher order contribution in which one of the jets is a gluon emitted from the quark 
box, that is 7(7 — > gq{q) or gq{q) with a spectator g or g of small px- We expect a more rapid 

fall-off from such a contribution. 

In order to calculate the Xg and distributions of Fig. 3 we integrate over the entire k"^ 
range of the gluon. The infrared contribution from k"^ < fcg is estimated using the prescrip- 
tion of ref. 1^. That is we use the strong-ordering approximation k'^ <C p\ and express the 
corresponding integrals in terms of the integrated gluon distribution g{xg, k^) at scale k^. This 
non-perturbative input is determined by the fit to the F2 data in a self-consistent way in the 
unified BFKL/DGLAP framework 0. 

In Table 1 we make the comparison of the data and the predictions of the dijet cross section 
for different values of the minimal pt of the jets. The calculation reproduces 70-80% of the 
observed rate. To put this comparison in context, we note that the Mepjet Monte Carlo predicts 
a dijet cross section for pT(jet) > 4 GeV of 2.8| or 2.6 nb according to whether GRV [|10] or 
MRS A partons were used [§. The latter set of partons have an integrated gluon more 
compatible with that used for our analysis. Thus the inclusion of small x contributions are 
seen to enhance the cross section, although the prediction of 3.2 nb is still below the measured 
value of 3.9 nb. 

^This number corresponds to Mepjet 2.0 with Q as the scale. If pT is taken to be the scale the cross section 
drops to 2.4 nb 1^. 
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r 1 


(J (expt) 


cr (theory) 


in GeV 


in nb 


in nb 


4 


3.9 


3.2 


5 


2.6 


1.8 


6 


1.6 


1.1 


7 


1.0 


0.7 


8 


0.6 


0.5 



Table 1: The comparison of the measured and theoretical integrated dijet production cross 
sections at the parton level for different pt cuts on the jet transverse momenta in the hadronic 
centre-of-mass frame. 

In principle dijet production appears to offer an opportunity to study the kx diffusion 
property of the BFKL gluon f{xg,k^). In practice the cuts, necessary on the transverse mo- 
mentum of the jets, curtail the small Xg "reach" of HERA, see (|^). Fortunately our unified 
BFKL/DGLAP approach is not restricted to small Xg, and gives a satisfactory description of 
the azimuthal decorrelations of the jets. Not surprisingly, the description is not unique and 
several standard Monte Carlos are known to also be able to accommodate the decorrelation 
data. The distribution, in the presently accessible kinematic domain, cannot therefore be 
regarded as a discriminator of the underlying small x dynamics. However our BFKL/DGLAP 
framework (with subleading ln(l/a;) contributions) does give an enhancement of the dijet rate, 
although the prediction still falls short of the observed cross section. Moreover by comparing 
the predicted dependence of the cross section with the data we are able to reveal the po- 
tential source of the remaining discrepancy. 

Acknowledgements 

We thank Brian Foster and Maciej Przybycien for discussions about the data. JK and AMS 
thank the Physics Department and Grey College of the University of Durham for their warm 
hospitality. AMS thanks Foundation for Polish Science for support. This work was supported 
in part by the UK Particle Physics and Astronomy Research Council, by the Polish State 
Committee for Scientific Research (grant no. 2 P03B 089 13) and by the EU Fourth Framework 
Programme TMR Network 'QCD and Particle Structure' (contract FMRX-CT98-0194, DG12- 
MIHT). 



4 



References 

[1] A.J. Askew, D. Graudenz, J. Kwiecinski and A.D. Martin, Phys. Lett. B338 (1994) 92. 
[2] J.R. Forshaw and R.G. Roberts, Phys. Lett. B335 (1994) 494. 

[3] S. Catani, M. Ciafaloni and F. Hautmann, Phys. Lett. B242 (1990) 97; Nucl. Phys. B366 
(1991) 135; S. Catani and F. Hautmann, Nucl. Phys. B427 (1994) 475; M. Ciafaloni and 



D. Colferai, [liep-ph/980635C 



[4] J.C. CoUins and R.K. Elhs, Nucl. Phys. B360 (1991) 3. 

[5] M. Przybycien, Ph.D. thesis "Two jets production in neutral current deep inelastic e'^p 
interactions at 300 GeV c.m.s. energy", Krakow, Dec. 1998. 

[6] V.S. Fadin, M.I. Kotskii and R. Fiore, Phys. Lett. B359 (1995) 181; V.S. Fadin, M.I. Kot- 
skii and L.N. Lipatov, [hep-ph/9704267| ; V.S. Fadin, R. Fiore, A. Flachi and M.I. Kot- 



skii, Phys. Lett. B422 (1998) 287; V.S. Fadin and L.N. Lipatov, |hep-ph/9802290 
Phys. Lett. B429 (1998) 127; V.S. Fadin, |heF"ph/9807527| , |tiep-ph/986752^ ; M. Ciafaloni 
and G. Camici, Phys. Lett. B386 (1996) 341; B412 (1997) 396; B417 (1998) 390 (E); 
Phys. Lett. B430 (1998) 349; M. Ciafaloni, |hep-ph/9709390[ 

[7] D.A. Ross, Phys. Lett. B431 (1998) 161; 
G.P. Salam, JHEP 9807 (1998) 019. 

[8] J. Kwiecinski, A.D. Martin and P.J. Sutton, Z. Phys. C71 (1996) 585. 

[9] J. Kwiecinski, A.D. Martin and A.M. Stasto, Phys. Rev. D56 (1997) 3991. 
[10] M. Gliick, E. Reya and A. Vogt, Z. Phys. C67 (1995) 433. 
[11] A.D. Martin, R.G. Roberts and W.J. Stirling, Phys. Rev. D50 (1994) 6734. 



5 



Figure Captions 



Fig. 1 Diagrammatic representation of dijet production in deep inelastic scattering via the 
photon-gluon subprocess. It defines the kinematics and dynamical quantities entering 
the factorization formula (^). 

Fig. 2 Theoretical predictions for the distribution with respect to the azimuthal angle between 
the q and q jets, compared to the data of ref. 0, for three different intervals of Bjorken 
X. The data are at the parton level and were obtained using the kr jet-finding algorithm. 
The predictions are normalised to the data point at = 2.55 radians. 

Fig. 3 Theoretical predictions for (a) the Xg dependence, and (b) the dependence, of the dijet 
production cross section compared to the parton- level data of ref. ^ . 
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